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The energy-zenith angular event distribution in a neutrino telescope provides a unique tool to 
determine at the same time the neutrino-nucleon cross section at extreme kinematical regions, and 
the high energy neutrino flux. By using a simple parametrization for fluxes and cross sections, we 
I present a sensitivity analysis for the case of a km'^ neutrino telescope. In particular, we consider 

, the specific case of an under-water Mediterranean telescope placed at the NEMO site, although 

' most of our results also apply to an under-ice detector such as IceCube. We determine the sensi- 

^S) ' tivity to departures from standard values of the cross sections above 1 PeV which can be probed 

independently from an a-priori knowledge of the normalization and energy dependence of the flux. 
We also stress that the capability to tag downgoing neutrino showers in the PeV range against the 
) cosmic ray induced background of penetrating muons appears to be a crucial requirement to derive 

^— ( , meaningful constraints on the cross section. 

PACS numbers: 13.15.+g, 95.55.Vj , 95.85.Ry 

P-t- I. INTRODUCTION 

O ^ 

^ I High energy neutrino astronomy is one of the most promising research lines in astroparticle physics. Neutrinos 
^ . are in fact one of the main components of the cosmic radiation in the high energy regime, and although their fluxes 
' are uncertain and depend on the specific production process, their detection would provide valuable information 
concerning the sources and the acceleration mechanism in extreme astrophysical environments. Similarly to photons 
and unlike cosmic rays, they keep directional information which can be used to perform astronomy. Differently from 
[ gamma rays, they are emitted only in hadronic processes and travel unimpeded through cosmological distances well 
1^ ' above the TeV energy. Moreover, due to the peculiar charge assignment of neutrinos in the electroweak standard 
model, they may be the particles most sensitive to exotic physics. An intriguing possibility is that a measurement of 
their cross section on ordinary matter in extreme and still unexplored kinematical regions might unveil new physics 
(see e.g. [H). 

From the experimental point of view, after the first generation of telescopes has proved the viability of the Cerenkov 
^ _ detection technique under deep water (Baikal 0) and ice (AMANDA 3) by detecting atmospheric neutrinos, one is 
, probably on the verge of the first detections at the IceCube Q telescope, being completed at the South Pole, and 
i;^ • possibly at the smaller ANTARES 0] telescope under construction in the Mediterranean. Additionally, ANTARES as well 
as NESTOR and NEMO Q are involved in R&D projects aimed at the construction of a km'^ Neutrino Telescope (NT) 
in the deep water of the Mediterranean sea, coordinated in the European network KMSNeT Q . 

In the present paper we address in a more quantitative way the question of how well a km"^ NT performs in the 
simultaneous determination of the neutrino flux and the neutrino-nucleon cross section by using as observable the 
energy-zenith angular event distribution. Note that, unless the spectrum is known a priori, for a NT surrounded 
by an isotropic medium the cross section and the flux would be completely degenerate. However, under the sole 
assumption that the incoming spectrum is isotropic, the different opacity of the Earth in different directions — which 
becomes relevant above ^100 TeV — allows one to disentangle the two observables, as it has been discussed before 
e.g. in [ol. [lol. . In other words, a NT is only sensitive to the product of cross section and flux, but the incoming 
isotropic flux is shielded (via a non-trivial function of the cross section) in a direction-dependent way, which makes 
possible to break the degeneracy by means of directional information. Here we study this property using a parametric 
and phenomenological approach, rather than illustrating it with a few discrete cases from models proposed in the 
literature. 

In particular, we consider the specific case of an under-water Mediterranean telescope, although most of our results 
also apply to an under-ice detector such as IceCube. Since the three different proposed sites for the under-water km'^ 
telescopes have shown event rate differences of the order of 20% [l^], for the sake of brevity we report the results of 
our analysis for the NEMO site only, which shows intermediate performances. 

This paper is structured as follows. In Sec. |TT]we introduce the parameterization used for the flux and cross section. 
In Sec. HIT] we describe the formalism used, while our results are reported in Sec. IIVI Finally, we conclude in Sec. |Vl 
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II. PARAMETERIZATION OF FLUX AND CROSS SECTIONS 



Many astrophysical objects are expected to produce a flux of high energy neutrinos, possibly with energies higher 
than the highest energy cosmic rays ever detected, of the order of 10^^ GeV. The center of mass energy available in 
a neutrino-nucleon collision for a neutrino with laboratory energy E^, is ~ \j2Ey m^r ~ 1.4TeVA/^^^/PeV, mjv 
being the nucleon mass. Already at PeV incident energies, this is well above the electro-weak symmetry scale of 
0(100) GcV and may thus reveal the onset of physics beyond the standard model^. Since neutrinos only respond to 
weak interactions, even weak effects may be revealed by a study of the neutrino-nucleon cross sections at NTs. One 
example is provided by low-scale quantum gravity models, see e.g. [T^ |. 

Unfortunately, there are two major problems to be addressed in such a program: the low statistics expected at 
NTs, which requires huge volumes and long operation times, and the ignorance of the diffuse astrophysical spectrum. 
At the range of energy Ey ~PeV which starts to be interesting for particle physics constraints, the natural target of 
NTs are diffuse fluxes rather than point-like sources. On the basis of existing predictions for extragalactic fluxes, we 
limit our analysis to energies above 10"°"'' PeV, where isotropic diffuse extragalactic fluxes are expected to dominate 
over the steeper atmospheric neutrino flux. Around PeV energies, the atmospheric neutrino event rates are very small 
even in a km'' telescope and the only remaining background is given by penetrating muons from extensive showers 
due to cosmic ray events, which may require an active veto to be disentangled from the signal. Given the observed 
behavior of cosmic ray and gamma ray spectra, it is reasonable to parameterize the neutrino spectra as a power law 
or, more generally, as a broken power-law. The latter parameterization would be probably needed when considering 
the flux over several decades of energy. However, for the following discussion most of the useful events will have 
primary energy Ei, < 10 PeV, and the single power-law approximation seems reasonable for a single decade in energy. 
The spectral index in Fcrmi-likc acceleration mechanisms is of the order of —2 or slightly steeper, but this observable 
is expected to be determined by the data as well. In the following, wc parameterize the neutrino flux (per flavor, 
summing neutrinos and anti-neutrinos) as 

where C and D are normalization and flux steepness free parameters, and for C = D = 1 one recovers the benchmark 
case of a Waxman-Balicall flux [l|, [l^ . 

Concerning the cross sections, we adopt a simple parameterization as a broken power-law where each energy interval 
(fixed a priori) has its own slope, and continuity is imposed at the boundaries. For example, the charged current 
(CC) cross section is given by 
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where we assume Ei = 10 "'^ PeV , E2 ~ I PcV and, for A = 1, the expression for the first energy bin reduces to the 
standard value reported in fT6| . An analogous parameterization is used for the neutral current (NC) cross section. 
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In the following, we shall consider both cases where CC and NC cross sections change proportionally to each other 
(i.e., A ~ A' and B = B') and the case where only ctnc is affected [A = B = 1 while A' and B' are set free). The 
former possibility mimics the scenario, for example, where the current extrapolation of parton distribution functions 
is wrong. It is more likely, though, that "truly exotic" new physics would manifest itself only in NC events; the 
latter case is a toy model for this class of situations. There are also cases where new physics manifests additionally in 
changes of the inelasticity of the cross sections, but we leave them outside the range of possibilities explored in this 
analysis: so, we shall make standard assumptions on the inelasticity of the collisions, as reported in more detail in 
the next Section. 



^ Of course, instruments looking at the highest energy events, say above 10^" GeV, have greater chances to detect significant departures 
from the standard model behavior, although the limited statistics likely requires future space based fluorescence detector (for an analysis 
of this case, see for example [Tsll ). 
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III. THE FORMALISM 

To perform our forecast, we follow the formalism outlined in [l^ . We generate neutrino tracks isotro pica lly outside 
the Earth and assume flavor equipartition due to neutrino oscillations and pion production dominance |17| . We place 
then a km^ NT at the NEMO site, see left panel of Fig. [T] For simplicity we assume the telescope fiducial volume as a 
cube bounded by six lateral surfaces, labeled D, U, S, N, W, and E for the Down, Up, South, North, West, and East 
side, respectively (see right panel of Fig. [T|). 

Despite the different behavior of the produced tau leptons with respect to muons in terms of energy loss and decay 
length, both i/^ and Vr event detection rates are sensitive to the matter distribution in the neighborhood of the NT 
site. In principle, the elevation profile of the Earth surface around the telescope may be relevant, similarly to what 
is known for Earth-skimming ultra- high energy i/^-'s at extensive air shower detectors, also investigated by some of 
the present authors [l^ for the Pierre Auger Observatory [l^ by using a Digital Elevation Map (DEM) of the site. 
In Ref. [13] the DEM's of the under-water Earth surface provided by the Global Relief Data survey (ET0P02) Q 
were used to estimate the effective aperture for Vr and i^f^ detection of a km'^ NT in the Mediterranean sea placed at 
any of the three locations proposed by the ANTARES, NEMO and NESTOR collaborations. It was found that the effect 
of the profile on the total number of events is negligible, although it may be important in the differential number of 
events from different directions, provided one has an excellent angular resolution. In the present analysis we take into 
account the underwater surface profile as done in Ref. p^ . although for the present results it plays a minor role. 



Moreover we also take into account the radial density profile of the Earth (reported, for example, in Ref. |21|). 

Neutrino telescopes were originally thought mostly as detectors, but they have some sensitivity to all neutrino 
flavors [l^l, although this does not trivially translate into a flavor tagging sensitivity. In particular, their capability 
as Vr detectors has become a hot topic in view of the fact that flavor neutrino oscillations lead to nearly equal 
astrophysical fluxes for the three neutrino flavors. While electron neutrinos only induce shower events inside the 
flducial volume, muon and very high energy tau leptons produced in charged current neutrino interactions create 
tracks which may be detected even when originating far from the instrumented volume. In the following, we consider 
as experimental observables: 

(i) The energy deposited in the detector, AE. 

(ii) The topology of the event, namely if it is a shower or a track event. 

(iii) Some (at least loose) information on the incoming direction of the event. 

In turn, these observables depend on some detector-dependent parameters like the optical properties of the medium, 
the geometry of the strings and the position of the optical modules on it, the efficiencies, etc. Here we take the 
simplified approach to consider a (muon or tau) track detected whenever the charged lepton decay length is longer 
than the length of the intersection of the trajectory with the instrumented cubic volume. Any other neutral current 
event which forms inside the instrumented volume, or charged current event which does not fulfill the previous 
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- LE TRACKS 

- - LE SHOWERS 

HE TRACKS 

- - - HE SHOWERS 




FIG. 2: In the left panel we give the zenith angle distribution of events of different kind and at different energy for the fiducial 
case A = B = C = D=1 and in one year of collection time. In the right panel we report the fractional contribution of the 
three flavors (empty, horizontal, and full hatching for electrons, muons, and taus, respectively) to the different type of events 
(T stands for tracks, S for showers), assuming flavor equipartition in the incoming fluxes. 



condition is classified as a shower. Although the actual criteria in a realistic experiment will differ, we expect this 
simplification to catch the physics of the problem. We shall comment on the impact of dropping the assumption of 
topological discrimination in Sec. IIVI 

Concerning the energy deposited in the detector, for our purposes it is sufficient to consider two bins only: a 
low-energy one (denoted LE) corresponding to 10"°'^ < A£'/PcV < 1 and a high-energy one (denoted HE) for 
A£^/PeV > 1. Note that we choose the lower bounds of the two bins equal to the quantities Ei and E2 in Eq.s ^ 
and Of course, this is a conventional choice dictated by simplicity and others are possible, although a posteriori 
it reveals to be a quite reasonable one. In order to calculate the energy released in the detector, one must also know 
the fraction of invisible energy and the energy loss mechanism of charged leptons. The fraction of energy not carried 
by the outgoing lepton and thus released to the struck nucleon is the inelasticity, y. At these energies y is basically 
independent of the flavor and is the same for i/ and P, decreasing from ~ 0.3 at E^, ~ i^i to ~ 0.2 at the highest 
energies of interest [2l|. If the interaction happens inside the detector (contained event) the energy yEi, is considered 
observed, since most of the energy in the hadronic shower generated by the struck hadron is released locally, virtually 
appearing as a point-like source^. For NC events, this is the only energy release, since the outgoing neutrino is 
invisible. For CC events, all of the remaining energy (1 — y) E^, is released locally for electron contained or other 
shower events, while for muon or tau tracks the energy losses are calculated as in [H, H^. Namely, the differential 
energy loss of the r leptons per unit of length in an underwater NT can be simply taken as di?^/dA = —(3rErQw, 
with = 0.71 X 1Q~^ cm^ g~^ and denoting the water density. Analogously, for muons one just needs to replace 
Pr with the corresponding value, = 0.58 x 10~^ cm^ g~-^. 

We end up with four different categories of events: LE tracks, LE showers, HE tracks, and HE showers. In the left 
panel of Fig. [2] we illustrate the zenith angle distribution of the events of different kind and at different energy for 
the fiducial case A = B = A' = B' = C = D~1 and in one year of collection time. First note that, in the case of an 
equal cross section in the two bins, the corresponding LE and HE curves should have the same shape, differing only in 
the normalization. The fact that this is not the case is a manifestation of the separate sensitivity of NTs to the cross 
section and the flux. The angular shapes of these curves are plausible, too. For downgoing events (cosi9 > 0) the 
interaction probability of neutrinos is ^1. Tracks can propagate into the detector even when forming relatively far 
from it; thus, the larger the grammage crossed the more numerous the events. The most notable manifestation of this 
effect is given by the peak near the horizon (cosi? ~ 0) due to the so called Earth-skimming events. The distribution 
of downgoing showers is instead flat. This is because contained events are the dominant contribution to this class, 
and the previous enhancement factor due to leptons created far away is not present. Instead, the distribution of 
the upgoing events (cosi? < 0) is more and more suppressed towards cos?? = —1. This is due to the fact that at 



^ Even the highest energy showers penetrate water or ice less than ~ 10 m, a distance short compared to any reasonable spacing of the 
photomultiplier tubes (PMTs). The radius over which PMT signals a re p roduced is ~250 m for a 1 PeV shower; this radius grows or 
decreases by approximately 50 m with every decade of shower energy [23| . 
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these energies the Earth is opaque to neutrinos, an effect most pronounced for events closer to the nadir. For tracks 
the suppression is even more pronounced, since it involves not only the suppression of the neutrino flux entering 
the detector, but also the additional energy loss for charged lepton tracks forming outside the instrumented volume. 
Comparing LE and HE events, note that HE events suffer a higher opacity, so the upgoing rate suppression is more 
pronounced with respect to the LE case. On the other hand, the range of charged leptons is larger, which explains 
the steeper behavior of the HE downgoing tracks curve. The number of downgoing track events also grows at HE, 
since the cross section is higher and we are in a regime where the interaction probability is ^ 1. The spiky bumps at 
cos'd ~ ±0.6 are simply an effect of the cubic geometry assumed (there is more instrumented volume when looking 
at angles closer to a vertex of the cube). In the right panel of Fig. [2l we show the flavor composition of the four 
classes of events. As expected, basically only muons contribute to the LE track bin. Indeed, the tau decay length 
is Dt ~ 50 m(£'T-/PeV). A realistic cut due to the spacing of the towers would probably remove even this small tau 
contamination. On the other side, the tau contribution to shower events decreases with energy, due to the increase of 
the relativistic tau gamma factor. Although we shall not discuss this point further here, it has been noted that the 
fact that different types of events are not flavor blind offers another possibility for astrophysical diagnostics [l^ . 

Although the overall number of events at a NT is limited, the previous discussion suggests that at very least a 
further partition of each of the previous bins in upgoing or downgoing events would greatly improve the diagnostic 
power for disentangling cross section from flux. Note that it is indeed a very weak requirement for a NT to be able 
to reconstruct at least the sign of the cosine of the zenith angle, in order to reject the background from downgoing 
muons from cosmic ray showers. Even so, it is at present unclear how well one can tag downgoing events as a signal 
against the cosmic-ray background. Around PeV energies, this should be feasible in IceCube, where the IceTop surface 
array offers a veto for cosmic ray events. For a Mediterranean km"^ neutrino telescope no final design exists, but the 
possibility has been suggested to extend the detection capability of the experiment with some sea top stations 
Keeping this possibility in mind, we shall consider both the case in which downgoing signal events can be used and 
the one where this is not possible. Also, we have tested that slightly changing the angular intervals and/or adding a 
third angular bin have only a minor effect on the sensitivity in the parameter space. 

For each choice of parameters, we end up with the counts iV/^, where K labels the type/topology of the event, 
i the energy bin and a the angular one. As already mentioned, we always consider two energy bins: a low-energy 
one, corresponding to Ei < AE < E2 and a high-energy one, corresponding to AE > E2- Regarding the choice of 
the angular binning, when both the downgoing and upgoing information is used, we divide the angular range in the 
two bins [0°,90°] and [90°, 180°] while, if only upgoing events arc included in the analysis, we consider the two bins 
[90°, 107°] and [107°, 180°]. We fix C = 1 since it is just a normalization, simply correlated to the exposure time 
needed to achieve the proper event statistics. The remaining parameters vary in a grid of 125-1-125 theoretical models 
as follows: D e [0.5, 1.5] in steps of 0.25, while A = A',B = B' e [0,4] in steps of 1 or A = B = 1,A',B' e [0,4] 
in steps of 1 in the two cases under investigation, that is: 1) CC and NC cross sections vary proportionally to each 
other; 2) only NC cross section changes. To determine approximately the range of parameters, we require that at 
least one event falls in each bin in one year of running. We then compare the counts N^^ with the expected counts 
C/^ for the benchmark case A = B = A' = B' = C = D = lhy means of a multi-Poisson likelihood analysis [1^ , in 
which the likelihood function, C cx cxp(— x^/2), is defined using the following expression for the x^: 



In this Section we present the results of our analysis. Fig.s [3] and |4] are obtained considering the case where CC 
and NC cross sections change proportionally to each other. The solid (black) curves in the panels of Fig. [3] show the 
marginalized constraints on couples of parameters (contours at 68% and 95% CL) in the case one has both upgoing 
and downgoing information available, for five years of exposure. Dashed (red) contours represent the case where 
no topological information is assumed, i.e. when one sums track and shower events. In general, there is a negative 
correlation between the cross section parameters A and B, since to some extent a higher cross section in the range 
El < E < E2 can be compensated by a smaller one in the higher energy range (note also that in the second line 
of Eq. ^ both A and B enter). Indeed, the primary energies are not directly observable , while the energy losses 
AE are. Also, there are negative correlations between A and C or D. The reason is that a higher cross section 
bin can be somewhat compensated by a smaller flux normalization (C) or a slightly steeper flux spectrum (D). The 
anti-correlation between A and B explains the correlations between B and C and B and D. In all cases the contours 
close, and the parameters can be determined to some accuracy. The degeneracy between flux and cross section is thus 
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FIG. 3: For the case where CC and NC cross sections change proportionally to each other, the solid (black) contours 
and 95% CL) represent marginalized constraints on pairs of parameters {^4, B, C, D} in the case where both upgoing and 
downgoing neutrino events can be tagged, for five years of exposure. Dashed (red) contours represent the case where no 
topological information is assumed, i.e. when one sums track and shower events. 



broken by the combined use of angular and energy independent information. Note how the topological information 
adds diagnostic power, although does not change the qualitative behavior. 

In Fig. [4] we show instead the marginalized constraints on couples of parameters (contours at 68% and 95% CL) in 
the case where only upgoing events can be tagged, again for five years of exposure. Solid (black) curves correspond 
to the case of two bins, [90°, 107°] and [107°, 180°], while dashed (red) curves show the results when these two bins 
are merged together. Only the flux parameters can be determined to some accuracy, even if some sensitivity to A and 
B is still recovered in the case of two angular bins. The basic physical reason for this behavior has been discussed 
analytically in [13]. As long as NC and CC cross sections change proportionally, then both the upward track rate 
and the upward shower rate are independent of cross sections, with or without new physics involved. If the energy 
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log 10 C 


0.04 ±0.10 


'-'•"^-0.19 


logio C 


-0 01+° " 


-0 01 +"■2" 

'-'■"^-0.24 


D 


1.0 ±0.1 




D 


1.1 ± 0.1 


1.1 ± 0.2 



TABLE I: Determination of cross section and flux parameters, obtained by marginalizing the (black) solid contours of Fig. |3] 
(left table) and of Fig. [5] (right table). 
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FIG. 4: For the case where CC and NC cross section change proportionally to each other, we report the marginalized constraints 
on couples of parameters {A, S,C, D} (contours at 68% and 95% CL) in the case where only upgoing neutrino events can be 
tagged and split in two angular bins, for five years of exposure. Track vs. shower separation capability is assumed. Dashed 
(red) contours represent the previous case where all events are summed in a single angular bin. 



dependence of the cross sections involved is the same, as we assume here for the fiducial case, the two observables 
(upgoing track and shower events) are just proportional to the integrated flux in the corresponding energy bin. 

The point we want to emphasize here is that without the tagging capabiUty of downgoing events, there is little 
hope for NTs to get meaningful constraints on the cross sections. A further binning of the upgoing events into two 
angular bins may only help a little bit, with the additional problem that the statistics is more limited in each bin. 

Finally, the curves in Fig.s [5] and [6] correspond to the solid (black) lines in Fig.s [3] and IH but assuming that only 
NC are changed (for illustrative purposes, we fix here A = B = \ for CC). The flux parameters C and D are quite 
accurately determined, in particular D is fixed at 10 % level, since the CC cross section contributing to tracks and, 
partially, to shower events is assumed to be known in this case. There is instead quite a poor determination of the 
NC cross section parameters, and in particular of A! . Indeed, only values of A! > 3.4 and B' > 2.4 are ruled out 
at 95% CL. This result is easy to understand. In fact, the largest contribution to NC shower events corresponds to 
downgoing neutrinos; so the event rate is almost independent of the intervening matter in the neighborhood of the 
NT site, being simply proportional to the product of the flux times cross section. For example, at high energies the 
number of NC events increases more than logarithmically as function of the highest possible neutrino energy: for 
B' > D/0.492 ~ 2 it becomes unacceptably large and thus excluded (actually the result would diverge if B' > 2 in 
the unphysical case of an incoming flux with D ~ I for arbitrary high energies.) The other relevant feature is the 
sensitivity only to a combination of A' and B' (roughly A' + kB', with k ~ 1.5 for the case of Fig. 5). The reason 
is that, due to the smaller energy transferred in visible channels in NC events, the NC events collected in both our 
energy bins are mostly sensitive to the high energy behavior of the cross section, i.e. the second line in Eq. (3). This 
function depends on a quasi-degenerate combination of the two parameters, explaining the results. Of course, a larger 
statistics at higher energies would eventually break this degeneracy, but would probably require too long exposure 
times. 

In Table[T]we summarize our results by reporting the parameter sensitivity at 68% and 95% after the marginalization 
of the likelihood function for both scenarios of CC and NC cross sections changing proportionally (left table) and 
only NC varying and standard behavior of CC (right table). The results are shown only for the best cases with both 
upgoing and downgoing tagged events and with track vs. shower events discrimination. 
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FIG. 5: For the case where only NC cross section changes while CC one assumes its standard value, we show the marginalized 
constraints on couples of parameters {A',B',C,D} (contours at 68% and 95% CL) in the case where both upgoing and 
downgoing neutrino events can be tagged, for five years of exposure. 




FIG. 6: Same of Fig. [5] but in the case where only upgoing neutrino events can be tagged and split in two angular bins, for 
five years of exposure. Track vs. shower separation capability is assumed. 
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V. DISCUSSION AND CONCLUSIONS 



111 this paper we have performed an analysis of the capabihty of a km"^ NT to disentangle the high energy diffuse 
neutrino flux and the neutrino-nucleon cross section in the PeV energy region for the incoming neutrinos. As well 
known, a separation is possible by exploiting the energy-zenith angular event distribution. Differently from previous 
treatments, we have used a simple phcnomenological paramctrization for fluxes and cross sections, presenting a 
sensitivity analysis in a given parameter space, rather than for a few benchmark cases. Our forecast is thus independent 
from strong assumptions on the normalization or energy dependence of the flux. The results we obtain are summarized 
in Fig.s m [21 m O ini and refer both to the cases where CC and NC change proportionally to each other, or the case 
where only NC events are affected by new physics. 

One important conclusion is that the capability to tag downgoing neutrino showers in the PeV energy range against 
the cosmic ray induced background of penetrating muons appears to be a crucial requirement to derive meaningful 
constraints on the cross section. In this case energy and zenith angle information on diffuse neutrino flux event rate in 
a NT (upgoing vs. downgoing) would greatly improve the diagnostic power for disentangling cross section from flux. 

Of course, our analysis presents several approximations. In particular, while we took into account the geometry of 
the site (for illustration we considered the NEMO site in the Mediterranean), we did not consider specific information on 
the detector geometry, efficiency, etc. Also, we assumed that in muon and tau tracks the energy loss is deterministic, 
rather than stochastic. A more accurate account of the telescope and of the energy losses is beyond the scope of 
this work. Yet, we expect that it would confirm our qualitative findings, while possibly worsening the accuracy with 
which parameters can be realistically reconstructed. Since our forecasts are likely more optimistic than a realistic 
treatment would show, our conclusion on the importance of downgoing event tagging is reinforced. In general, the 
background due to single energetic muons generated from cosmic rays and surviving deep underground drops faster 
with energy than the expected neutrino signal, and well above the PeV region the contamination should be negligible. 
However, detailed simulation are required to settle down this issue quantitatively; an active veto may help pushing 
down in energy the region where the background can be disentangled from the signal on an event-by-event basis. 
For the IceCube telescope, in the PeV energy range downgoing muons from cosmic ray showers should be vetoed 
thanks to the IceTop array. For the Mediterranean km'^ neutrino telescopes no final design is available, yet. Given 
the importance of the present topic, we believe this is an important issue to be addressed in the design phase of these 
telescopes. 

Even in the case where no downgoing tagging for neutrinos will be available for a Mediterranean NT, still its upgoing 
events at the PeV would provide important constraints on the flux and cross section if a combined analysis of its data 
with IceCube ones is performed. At very least, it should allow to test the hypothesis that the diffuse neutrino fluxes 
at the PeV from the two hemispheres are consistent with each other, thus validating the hypothesis of an overall 
isotropic flux. 
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